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We consider the hypothesis that dark energy and dark matter are the two faces of a single dark 
component, a unified dark matter (UDM) that we assume can be modeled by the afhne equation 
of state (EoS) P = po + op, resulting in an effective cosmological constant pA = — po/(l + o)- The 
affine EoS arises from the simple assumption that the speed of sound is constant; it may be seen as 
an approximation to an unknown barotropic EoS P — P{p), and may as well represent the tracking 
solution for the dynamics of a scalar field with appropriate potential. Furthermore, in principle 
the affine EoS allows the UDM to be phantom. We constrain the parameters of the model, a and 
f^A, using data from a suite of different cosmological observations, and perform a comparison with 
the standard ACDM model, containing both cold dark matter and a cosmological constant. First 
considering a flat cosmology, we find that the UDM model with affine EoS fits the joint observations 
very weU, better than ACDM, with best fit values a = 0.01 ± 0.02 and Ha = 0.70 ± 0.04 (95% 
confidence intervals). The standard model (best fit SIa ~ 0.71 ± 0.04), having one less parameter, 
is preferred by a Bayesian model comparison. Ifowever, the affine EoS is at least as good as the 
standard model if a flat curvature is not assumed as a prior for ACDM. For the latter, the best fit 
values are Qk = -0.02j;°:°2 

and JIa ~ 0.71 ± 0.04, i.e. a closed model is preferred. A phantom 
UDM with affine EoS is ruled out well beyond 3a. 

PACS numbers: 98.80.-k; 98.80.Es; 95.35.+d; 95.36.+X 



I. INTRODUCTION 

In the past few years, evidence that the universe is 
going through a phase of accelerated expansion has be- 
come compelling [1-4]. The cause for the acceleration, 
however, remains mysterious, and might be regarded as 
the most outstanding problem in contemporary cosmol- 
ogy. The favoured working hypothesis is to consider a 
dynamical, almost homogeneous component with nega- 
tive pressure, dubbed dark energy [5-8]. Such a frame- 
work helps alleviating a number of fundamental problems 
arising when an ad hoc cosmological constant A term in 
Einstein equations is interpreted as the energy density of 
the vacuum [9-11]. Alternatives include modifications of 
gravity and extra dimensions [12-14]. 

Here, within the framework of General Relativity 
(GR), we explore the hypotesis that dark energy (the 
source for the observed acceleration of the universe) and 
dark matter (DM, required to explain structure forma- 
tion) are the two faces of a single dark component, a 
unified dark matter (UDM) that we assume can be mod- 
eled by a simple 2-parameter barotropic equation of state 
(EoS). In turn, this EoS can be derived from the simple 
assumption that the UDM speed of sound is constant. 

Recent analyses of cosmological data suggest that 
within the standard ACDM scenario (Cold DM plus A) a 
slightly closed, positively curved model is preferred [15]. 
With this in mind, motivated by a theoretical bias in fa- 
vor of a flat inflationary scenario, we compare our UDM 
model, assumed to be flat, with both flat and curved 



ACDM models, using Bayesian methods [16, 17]. 



II. THE MODEL 

Let's consider a flat Friedmann-Robertson- Walker uni- 
verse in GR, with radiation, baryons and a single UDM 
component with energy density px ■ The dynamics of this 
model is governed by the Friedmann equation 

H^ = —{Pr+Pb+Px) (1) 

and the energy conservation equations of the three com- 
ponents. Here H is the Hubble expansion scalar related 
to the scale factor ahy H = d/a, and no ad hoc cosmo- 
logical constant A term is assumed in Eq. (1). 

With the usual scaling laws pb oc and pr cx: a^'^ 
for the density of baryons and radiation, we now assume 
that the dark component is represented by a barotropic 
fluid with EoS Px = Px{px), satisfying the conservation 
equation 

PX = -3H{px + Px). (2) 

It is clear from Eq. (2) that, if there exists an energy 
density value px = Pa such that Px{pa) = — PA, then 
PA has the dynamical role of an effective cosmological 
constant: pa — (see [18] for a detailed discussion). 

In order to provide for acceleration, our UDM com- 
ponent must violate the strong energy condition (SEC) 
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(see e.g. [19]): Px < — Px/S at least below some red- 
shift. This can be achieved by a constant wx = Px/pxi 
which is indeed allowed by observational tests based only 
on the homogeneous isotropic background evolution, see 
e.g. [20, 21] (tests of the same kind we are going to con- 
sider here, see below). In this case, however, one would 
have Cg = dPx/dpx = w < for the adiabatic speed of 
sound, and this would have nasty consequences for struc- 
ture formation in an adiabatic fluid scenario"'^. Instead, 
given that the EoS Px = Px{px) is unknown, the next 
simplest approximation we can make to model it is to as- 
sume a constant speed of sound dPx/dpx — a, leading 
to the 2-parameter affine form [18] 

Px'^Po + apx- (3) 

This allows for violation of SEC even with = a > 0. 
Then, using (3) in (2) and asking for pa — leads to the 
effective cosmological constant pA = — Po/(l+a)- Eq. (3) 
may also be regarded (after regrouping of terms) as the 
Taylor expansion, up to 0(2), of any EoS Px = Px{px) 
about the present energy density value pxo [25] . In addi- 
tion, it may also represent the exact tracking solution for 
the dynamics of a scalar field with appropriate potential 
[26]. 

The EoS (3), if taken as an approximation, could be 
used to parametrize a dark component (either UDM or 
DE) at low and intermediate rcdshift. In the following 
we are going to make a more radical assumption, that 
is we are going to extrapolate the validity of Eq. (3) to 
any time. In doing this we are therefore going to build 
a cosmological model based on a UDM component with 
EoS (3), to be tested against obscrvablcs at low and high 
redshift, as described in the next section. 

Using the EoS (3) in the conservation equation (2) 
leads to a simple evolution for px{o)'- 

Pxia) = PA + {pxo-PA)a-^^^+"\ (4) 

where today Co = 1. Formally, with the EoS (3) we can 
then interpret our UDM as made up of the effective cos- 
mological constant pA and an evolving part with present 
"density" pm = pXo — Pa- We may then dub it AaDM^. 
A priori, no restriction on the values of a and po is re- 
quired, but one needs Po < and a > — 1 in order to 



^ The EoS with constant w may as well represent a scalar field 
with exponential potential [22], and this would have a unitary 
effective speed of sound, c^ff = 1, see e.g. [23, 24] and references 
therein. 

^ It should be clear that, for pm > and starting from = a = 
const., our model is formally totally equivalent to a model where 
a standard cosmological constant A = SnGpA is assumed a-priori 
in Einstein equations, i.e. in (1)), together with a standard fiuid 
with linear EoS w = a and present density pm- Indeed, by 
definition pjy = 0, thus it is easy to check that a-nd pm = px ~ 
PA in (4) separately satisfy their own conservation equations, so 
that there is no coupling between the dark energy part pA a^nd 
the dark matter part pm of our UDM fluid. 



satisfy the conditions that pA > and p ^ pA i'a the 
future, i.e. to have that pA is an attractor for Eq. (3). 
Then, our UDM is phantom if pm < 0, but without a "big 
rip", cf. [18]. Last, but not least, it follows from Eq. (4) 
that, with pm > and a = 0, our model is equivalent to 
a ACDM. This allows for a straightforward comparison 
of models, done in section IV. 

Finally, the EoS of our AaDM model can be charac- 
terized by its parameter wx = Px/ Px'- 

wx = -{l + a)^+a, (5) 

PX 

wxo - -nA + a{i-nA}-nA{i-nxo) + --- (6) 

where today's value wxo is approximated using flxo ~ 1 
and assuming |a| <^ 1. Hence at leading order wxo — 
— Oa, while Wx — > —1 asymptotically in the future. 

III. OBSERVABLES 

In order to constrain the affine EoS (3) we perform 
maximum likelihood tests on its parameters, using differ- 
ent cosmological probes. We assume a flat cosmological 
model, n ~ 1, and we hold the baryon density fixed to the 
best fit value derived from the analysis of WMAP 3 year 
data fib/i^ = 0.02229 ±0.00075 [15]. We also fix the Hub- 
ble parameter to the value measured by the Hubble Space 
Telescope (HST) Key Project Hq = 72±8 km/s/Mpc [27] 
(but see below about marginalization). This leaves only 
two free parameters in our analysis, i.e. those character- 
izing the affine EoS: the constant speed of sound c^ = a 
and flA = SttGpa/ (3Ho)- Given that our model contains 
an effective cosmological constant parametrized by Qa 
and that for > in (4) and = a = this model is 
equivalent to a ACDM, in essence we are going to con- 
train with observables the possible variations of the speed 
of sound of the overall dark component from the value it 
takes in the ACDM case, under the simple assumption 
that Cg = a is a constant. 

The observables we use have become the standard tools 
to probe the background cosmology, through their de- 
pendence on the expansion rate of the universe (see, e.g. 
[20, 21]). They are: 

a. The present age of the universe, to- This is ob- 
tained directly from the definition of H. We compared 
the theoretical prediction on to with its best estimate 
to = 12.6t.2A Gyr, derived from a combination of differ- 
ent astrophysical probes in [28] . 

b. The luminosity distance of type la supernovae. 
Type la SNe light curves allow a determination of an 
extinction-corrected distance moduli, 

p.o = m- M = 5\og{dL/Mpc) + 25 (7) 

where dL = {L/A-KP f/'^ = [l + z) J^dz'/H{z') is the 
luminosity distance. We compare our theoretical predic- 
tions to the values of fiQ estimated in [4] using a sample 
of 182 type la SNe (the new Gold dataset [4] containing 
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the old Gold sample [3], the Supernovae Legacy Survey 
(SNLS) [29] and 12 new SNe observed by the HST). 

c. The location of cosmic microwave background 
(CMB) acoustic peaks. When the geometry of the uni- 
verse is held fixed at = 1, this only depends on the 
amount of dark matter, through the shift parameter TZ 
defined as: 



n = n]l^HoDA{zis), Da{zis) = 



dz' 



(8) 



where Da is the distance to the last scattering sur- 
face. We have identified pm = Pm for our model, i.e. 
fim = 87rG'/5m/(37Jo), and the redshift of last scatter- 
ing has been estimated using the fit function from [30], 
which for our choice of fift/i^ gives zis = 1089. We derived 
a value of TZ using five Monte Carlo Markov chains pro- 
duced in the most recent analysis of the WMAP data^. 
We estimated a value TZ = 1.71 ±0.03. Our new analysis 
updates previously pubhshed values of 7?., e.g. [31]. 

d. Baryon acoustic oscillations. The recent detec- 
tion of acoustic features in the matter power spectrum 
measured with the Sloan Digital Sky Survey (SDSS) 
galaxy survey constrains the A parameter defined as: 



A 



Ho 

0.35 c 



Dv = 



D\{z) 



H{z) 



(9) 



where Dy is the distance to z = 0.35 taking into account 
the distortion along the line of sight due to the redshift. 
The value of the A parameter adopted in our analysis 
is the one measured from the SDSS luminous red galaxy 
survey: A = 0.469 ± 0.017 (for = 0.98) [32]. 

We perform our likelihood analysis by sampling the 
parameter space in the range [—0.3, 0.3] for a (where the 
upper limit is due to the requirement that the UDM fluid 
does not scale faster than radiation) and [—0.3, 0.9] for 
r^A- The size of our parameter space does not pose any 
serious problem in terms of computational time. Thus, 
we preferred to perform a direct evaluation of the likeli- 
hood on a predefined grid in the volume under considera- 
tion, rather than adopting the now popular Monte Carlo 
Markov chain approach. When calculating confidence in- 
terval in the a — i^A plane from the SNe data, we always 
marginalize over a calibration uncertainty (treated as a 
nuisance parameter) on the absolute magnitude. Note 
that marginalizing over calibration uncertainty is com- 
pletely equivalent to marginalizing over TJqj since both 
parameters effectively act as an additional term in the 
distance modulus (7). 

A final caveat on the use of the R and A parameters. 
The numerical values adopted in our analysis were both 
derived assuming an underlying ACDM cosmology. In 
principle, one might question whether it is a good as- 
sumption to adopt the same numerical values when con- 
straining a different class of models. We believe that this 



is a reasonable approximation when comparing the evo- 
lution of the expansion rate with observations, and one 
that has been adopted in many other previous analyses 
of non-standard dark energy models [20, 21, 33, 34]. 

Indeed, the physical motivation for adopting this ap- 
proximation lies in the fact that both observables mostly 
depend only on the expansion rate of the universe. The 
CMB peak position is defined by a standard ruler (the 
sound horizon at recombination) which is very weakly 
affected by the dark energy component (firstly, because 
dark energy has essentially no influence before recombi- 
nation, secondly, because the sound speed which enters 
into consideration here is that of the baryon-photon fluid, 
regardless of the other components). This standard ruler 
is then seen under different angles depending on the an- 
gular diameter distance at recombination. This, in turn, 
depends only weakly on the equation of state of dark en- 
ergy, and more strongly on its density: in our model, at 
least for the values of the a parameters considered here, 
this follows very closely that of ACDM, hence it is rea- 
sonable to believe that the approximation does not intro- 
duce significant errors. The same considerations apply to 
the A parameter. On this basis, we are then confident 
that the errors we are making by extending the ACDM 
approximation to our model lie well within the present 
accuracy of our analysis. 



IV. RESULTS 




FIG. 1: UDM likelihood contours at 68%, 95% and 99% C.L. 
(shaded regions), in the a-Q,A plane, obtained from the age 
of the universe, type la supernovae (SN la), the location of 
cosmic microwave background acoustic peaks (CMB), and 
baryon acoustic oscillations (BAO). Also shown (continuous 
curves) are the likelihood contours obtained from combining 
the datasets. 



^ See the Lambda web site: http://lambda.gsfc.nasa.gov. 



Fig. 1 shows the confidence levels derived from the dif- 
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FIG. 2: Likelihood contours at 68%, 95% and 99% C.L. in 
the a — f^A plane, obtained from the combined datasets. The 
dot represents the combination of parameters which best fits 
our data. The vertical line locates the models describing flat 
ACDM models. 



ferent cosmological observations used in our analysis, to- 
gether with the joint confidence levels from the combined 
datasets. The combined levels are also shown separately 
in more detail in Fig. 2. A fairly narrow area of the pa- 
rameter space is identified by our analysis, encompass- 
ing models which are close to the standard cosmological 
constant case with Ha ~ 0.7. However, there is a pref- 
erence for values of a larger than zero. Our best fit is 
= 0.70 ± 0.04 and a = 0.01 ± 0.02, both at 95% 
C.L.; for these values, we find = 157.64 with 185 data 
points. The marginalized likelihood functions for the 17a 
parameter are shown in Fig. 3. As said above, when 
a = is assumed as a prior, our model describes a flat 
ACDM model. In this case, we find a result which is con- 
sistent with previous analysis: flA = 0.71 ± 0.04 at 95% 
C.L., with = 158.83. We note that, for a = 0, the 
SNe data seem to prefer a slightly smaller value of fl\ 
(see Fig. 3, top panel). This agrees with similar results 
(see e.g. [35], which also suggest that this might be due 
to some systematic effect in the latest SNe data). 

In order to assess the statistical significance of our find- 
ings, we compare the AaDM model with the fiat ACDM 
using the Akaike and the Bayesian Information Criteria 
(AIC and BIC, [16]), as well as the Bayesian evidence E 
[17] (see also [36, 37]). In essence, all these can be used 
as model selection criteria, i.e. they estimate how much 
adding a new parameter to the cosmological model (in 
our case, the a parameter) is justified by the increased 
goodness of fit (cf. also [20]). BIC improves on AIC tak- 
ing also into account the number of data points. BIC 
should be a crude approximation to where the latter 




FIG. 3: Marginalized likelihood functions for Qa derived from 
each dataset and from their combination. Top panel: the flat 
ACDM model, obtained when q = is assumed as a prior. 
Bottom panel: the likelihood obtained by assuming the best 
fit value a = 0.01 for the a parameter. 



is a further refinement as it averages the likelihood of the 
model in the prior: E = J C(j>)P{p)dp where p are the 
free parameters of the model and P is the prior distribu- 
tion of the parameters. E is precisely the likelihood of 
the model given the data, according to Bayes theorem, 
and as such gives positive weight to a "good on average" 
fit over a larger volume in parameter space. The model 
comparison is achieved through computing the quantity 
A\nE = IuEao - Ih^^aq- We find AlnE ^ 3.9, which 
has to be regarded as a rather strong evidence in favor of 
the fiat ACDM model, given the present data [17]. Us- 
ing AIC and BIC gives milder results (evidence but not 
strong evidence, cf. [16]), summarised in Table I. 

However, our analysis shows rather clearly that the 
best fit model in our class is not a fiat ACDM and, as said 
in the introduction, recent results favor a slightly closed 
universe within the ACDM models [15]. Then, a question 
worth investigating is which model performs better in 
terms of AIC, BIC and E if we add an extra parameter 
to the flat ACDM, either by allowing for curvature in 
the ACDM scenario, or through our AaDM models, i.e. 
adding or a, respectively. We find that for non-fiat 
ACDM models the best fit is achieved for JIa = 0.71 ± 



0.04 and 



-0.02 



+0.01 
-0.02' 



see Fig. 4, and again there 



is evidence in favor of the flat ACDM, see Table I. A 
direct comparison of the not nested [16] non-fiat ACDM 
and the AaDM models is possible in terms of AIC and 
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BIC: in terms of the latter, we see from Table I that 
there is positive evidence in favor of the AaDM model. 
An indirect comparison of the AlniJ's suggests however 
that the two models are statistically almost equivalent. 



Model 




AAIC 


ABIC 


AlnB 


flat ACDM 


158.83 











flat AqDM 


157.64 


0.8 


4 


3.9 


curved ACDM 


161.53 


4.7 


7.9 


3.4 



TABLE I: Model comparison with information criteria and 
Bayesian evidence: the A's compare the flat AaDM and the 
non-flat ACDM against the standard flat ACDM model. 
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FIG. 4: ACDM models: likelihood contours at 68%, 95% and 
99% C.L. in the fix — JIa plane for the combined dataset. The 
dot represents the combination of parameters which best fits 
our data. The vertical line locates the flat (fix = 0) models. 



V. CONCLUSIONS 

Today's standard model for the homogeneous isotropic 
universe is the 1-parameter flat ACDM, based on GR as 
the theory of gravity, a flat spatial geometry as it fol- 
lows from the inflationary paradigm, cold dark matter 
and an ad hoc cosmological constant A as the only two 
relevant constituents apart from baryons and radiation. 
However, the now well established accelerated expansion 
of the universe calls for an explanation beyond the ad 
hoc A term in Einstein equations. In the game of mod- 
els testing and selection, in absence of really compelling 
models from fundamental physics, that of simplicity of 
models - the Occam's razor - remains a critical ingredient. 
Quantitatively, this can be implemented using Bayesian 



methods in the statistical analysis of data, thereby mea- 
suring the worth of models by favoring those that give 
a good fit with fewer parameters. On the other hand, 
it does makes sense in our view to consider models built 
on some physical basis, and compare models with equal 
number of parameters but different physical motivations, 
in an attempt to establish which model performs better. 

Here we have assumed a single dark component with 
barotropic EoS Px — Px{px) to make up for a total 

= 1 (i.e. a flat model) and the acceleration of the 
universe. Assuming a constant speed of sound c^ — a the 
2-parameters affine EoS (3) then follows: we then dub our 
model as AaDM, as it admits an effective cosmological 
constant and thus it is formally equivalent to having a 
A term and an evolving part characterized by a. For 
a = 0, our UDM is formally equivalent to ACDM for 
today's density value pxo > PA in (4). On the other 
hand, for any a > — 1 and for pxo < PA our UDM would 
be phantom, with no "big rip" [18]. 

Simple parameterizations of the EoS may lead to mis- 
leading results [4, 38]. Comparing 2-parameter EoS's, 
in our view barotropic forms^ like (3) have at least the 
advantage over the common scale factor or redshift pa- 
rameterizations of relating directly Px and px- This 
seems more physical for an EoS: one wants the param- 
eterization to express an intrinsic and realistic relation 
between the fluid or field [39] variables, disentangling it 
from universe evolution. In addition, an EoS of the form 
Px ~ Px[px) has a vahdity that goes beyond that of 
the homogeneous isotropic background. It can indeed be 
easily carried forward to study inhomogeneities: at per- 
turbative order one can use it assuming adiabatic pertur- 
bations with 5p — c^Sp, representing either a barotropic 
fluid or a k-essence field [26] , or else to formally describe a 
scalar field with unitary effective speed of sound, c^y^ = 1 
[23, 24, 26]. Wc are addressing the growth of perturba- 
tions and their CMB signatures for all these cases in a 
forthcoming paper [40] . 

Constraining the AaDM model with various observ- 
ables we find that it fits the data well, better than 
ACDM, with best fit values a = 0.01 ± 0.02 and Qa = 
0.70 ± 0.04. From Eq. (6) we then get wxo - -0.7 at 
leading order: indeed we find that values of a and Q\ 
corresponding to a phantom dark component are ruled 
out well beyond 3cr's. Our results for model comparison 
are summarized in Table I. The flat ACDM model is still 
preferred, with a "rather strong" evidence [17]. Com- 
paring our model with a non-fiat ACDM shows however 
that the two are, all in all, almost statistically equivalent. 
Then, one would advocate our model, which is assumed 
to be flat, on the basis of an inflation motivated theoret- 
ical bias. 

Beyond this factual analysis, we may as well look at 



* See e.g. [20] for a recent observational test of the Chaplygin gas 
and other examples. 
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Fig. 1 from a different perspective, it is indeed rather 
clear that BAO and SNe data alone prefer a larger = a 
value than our best fit value. Indeed, we find that the 
BAO+SNe best fit is a « 0.1. On the other hand, it is 
also clear from Fig. 1 that the CMB strongly constrains 
a, and forces it to take a very small value (cf. [40, 41]). 
This seems to suggest that high redshift data (CMB) 
favor Cg ~ 0, while low and intermediate redshift data 
(BAO and SNe) allow for a larger c^, i.e. that the speed 
of sound of a UDM component should be variable. This 
is indeed in line with the findings of an analysis of models 
such as the generalized Chaplygin gas, which shows in- 
deed a very good fit to data [20]. As we have shown here, 
the standard ACDM can be seen as a UDM with vanish- 
ing speed of sound at all redshifts, = a = 0. Thus, If 
proved by future observations, a non vanishing speed of 
sound at low and intermediate redshift for the whole dark 
component in the Friedmann equation (1) would provide 
a sort of observational "no-go theorem" for the standard 
ACDM in GR. 

In this paper, in order to consider a minimalist 2- 
parameter model, we have assumed a single UDM compo- 
nent with constant speed of sound = a in addition to 
baryons and radiation. It could well be that the absence 
of coUisionless CDM is going to spoil structure formation 
in this model, or else, that our model is subject to much 



stronger constraints, forcing it to coincide with ACDM 
in practice, when a full CMB analysis is carried out (cf. 
[40, 41]). On the other hand the simple AaDM compo- 
nent, derived assuming a constant speed of sound, may 
turn out to be a useful and physically motivated model 
for dark energy (cf. [42, 43]), additional to CDM. In any 
case, it can be shown [26] that the affine EoS (3) cor- 
responds to an exact solution for a quintessence scalar 
field with appropriate potential, as well as to the general 
dynamics of a k-essence field. The quintessence field in 
particular is going to have different perturbations than 
the corresponding fluid (see e.g. [23, 24]), which then 
need to be separately analized. All these problems can 
only be settled with future work. 
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